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Abstract: The CP violating asymmetry between the partial decay rates of W^h? 
and W~hP is calculated in first order of the weak coupling constant a^^ = /^ir in 

MSSM with complex parameters. The dependence on the phases of At and Mi is discussed. 
Different values of tan/3 are considered. The asymmetry is up to the order of 10~^. 
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1. Introduction 



The Minimal Supersymmetric Standard Model (MSSM) implies the existence of a pair of 
charged Higgs bosons H^. At tree level there are three possible decay modes of into 
ordinary particles: tb, ut^ and W^hP where is the lightest neutral 

Higgs boson. The lighter fcrmions from the first two generations and the heavier neutral 
Higgs bosons are not considered. Loop corrections due to a Lagrangian with CP violating 
phases lead to decay rate asymmetries between the partial decay widths of and H~ and 
that would be a clear signal of CP violation. In refs.[l] and [2] such decay rate asymmetries 
are considered in the MSSM with complex phases for the quark decay mode H'^ tb - the 
asymmetry 5tb, and for the Icpton decay mode ^ vt - the asymmetry 8ut- In order 
to finalize this investigation we consider here the decay rate asymmetry of —>■ W^hP: 

_ r{H+ W+h'^) - r{H- w-h'^) 

~ r{H+ ^ w+ho) + r{H- ^ w-W) ■ ^ ^ ^ 

We shall work in MSSM with complex parameters in first order of the weak coupling 
constant a^, = 5^/47r. Within MSSM, after redefining the fields, the new sources of CP 
violation are the phase of the higgsino mass parameter n, two of the phases of the gaugino 
masses Mj, i = 1,2,3 (we choose these to be the phases of Mi and M3), and the phases 
of the trilinear couplings of the fermions f,Af. Especially the latter ones are practically 
unconstrained. Previously this asymmetry was considered in the two-Higgs doublet model 
in [3]. 

Discussing CP violation in the decay widths, we must keep in mind the branching ratios 
of the relevant decay modes. vt~^ is significant for low mjj+, below the tb threshold. 

This determines the sensitivity of Si^r to light i>f and x^x^ in the loops, i.e. to the phases 4>t 
and (f>i of At- and Mi . The decay — > tb dominates for high mH+ , which determines the 
sensitivity of Stb to the phases of At and Ai,. The complication with the decay H'^ — W^h^ 
is that the final state hP is not observed yet and m^jO is an unknown parameter. However, 
once and tan /3 are fixed, the SUSY structure of the theory determines uniquely both 
m/jO and the coupling H^W~h^. There are two consequences that are important for us. 
First, increasing mff+, the mass m/jO is saturated approaching its maximum value. At 
tree level this is particularly simple: m/jO < m^"^ = nizl cos2/3| [4], while including QCD 
and SUSY radiative corrections m^"^ can be increased considerably, mj^"^ c± 130 GeV 
(for a recent review see [5] and the refs. therein). There is also an experimental lower 
bound for m^o, m^o > 96 GeV [6]. Thus, respecting both the experimental and theoretical 
bounds, we shall consider m^o in the range 96 < m^o < 130 GeV. In this range of mj^o, 
for 171^+ > rnw + mj^o we are already in the saturation regime where we may keep 
and m^o as independent parameters. The second consequence concerns the H'^W'h'^ 
coupling which determines the Br(ii"+ — W^h^). Increasing it quickly falls down 

and, depending on tan/?, we can enter the so called decoupling limit, cos^(/9 — a) — ^ 0, 
where the Br(i7+ W^h^) almost vanishes. This imposes severe restrictions on mfj+ 
and tan (3. In order to keep the value of Br{H~^ W~^h^) at the level of few percents, we 
shall consider 200 < mj^+ < 600 GeV and low tan/3, 3 < tan/? < 9 (tan/? < 3 is already 
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excluded from the Higgs searches at LEP[6]). In accordance with this, will receive 

contributions from i>f^, tb and x^x^ in the loops. Due to the large mass of the top-quark, 
the radiative corrections with stops and sbottoms with low masses appear to be too large 
to be considered within the approximation used here. This means that we assume that 
the squarks are heavy and will not contribute in the considered range of . This will be 
discussed in the next Section 2. Thus, we shall consider the sensitivity of the asymmetry 
6wh to the phases 0r and (pi of the chargino-neutralino and the slepton sectors. According 
to the experimental limits on the electric dipole moments of the electron and the neutron, 
we assume a zero phase for the Higgsino mass parameter n, (p^ = 0. However, as is shown 
in [7], a large phase is not impossible, it would require fine-tuning between the phase 
(pn and the other SUSY parameters. We end up with a short discussion of the influence of 
(j)^ on 6who. 

The paper is organized as follows. In the next section we present the analytic expression 
for the asymmetry. In Section 3 we discuss the numerical results in MSSM. We end up 
with a conclusion. There are two Appendicies - with the Lagrangian and with the analytic 
expressions for the imaginary parts of the Passarino Veltmann (PV) integrals. 



2. The asymmetry 

We write the matrix elements of H'^ — ^ W^hP and H" — ^ W~hP in the form: 

MH±=ige^{pw)pty^, (2.1) 
where e^{pw) is the polarization vector of W"^, Y"^ are the loop corrected couplings: 

Y^ = y + 6Y^ + 6Y^ + 6Y^ + .... (2.2) 
Here y is the tree level coupling: 

y = cos(a — /3), (2.3) 

and SY^,k = 1,2,3,... are the SUSY-induced loop corrections. 
The decay rates of — > W^hP are: 



lb 'mH''^w 



r{H^ W^h") = v-^' -/»'-i^^ |y±|^ (2.4) 

Here: 



A(x, y, z) = + y"^ + — 2xy — 2xz — 2yz, (2-5) 

|y±|2 =y2 + 2y^5ie(5y±)+0(Q^). (2.6) 
k 

At tree level the decay widths of and H~ are equal and there is no CP violation: 



r,ree{H^ W^h') = ^'"^ ' ^""^^ cos^a - (3). (2.7) 

16 m^H'^w 



A3/2(rn2^, 77? 2, rn2) 
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CP violation is induced by loop corrections. They have CP-invariant and CP-violating 
contributions: 

5Y^ = SYl""" ± 5Y^^ (2.8) 

and each one has real and imaginary parts: 

^yinv ^ sfie((^yfe*™) + i^mCc^yfe^""), SY^^ = ^e{5Y^P) + &m{5Y^P) (2.9) 

The asymmetry is determined by '^e{5Y^^). pTirthcr we shall work in first order 
of the weak radiative coupling constant = /Att. This approximation means that we 
neglect the CP-invariant loop corrections 'ike{5Y^^'") in the denominator in formula (1-1). 
Then from (1.1), (2.4) and (2.6) we obtain: 

6who ^ — . (2.10) 

Here the sum is over the loops with CP-violation. We work in MSSM with complex 
parameters. The Lagrangian is in Appendix A. 

In general, there are two types of SUSY radiative corrections - self-energy loops and 
vertex corrections. We have self-energy loops on the VF+-line, on the hP-linc and on the 
i7+-line. The loops on the VF^-line are proportional to the 4- momentum of W^, py^, and 
because of the gauge invariance {p^^eaipw) = 0); equal to zero. As S-^y^o is determined 
by the absoptive parts of the loops, the loops on the /i'^-line will contribute only if the 
kinematic condition m^o > ("^-i + "^2)^ is satisfied, where fhi and 777-2 are the masses of 
the two particles in the loop. However, because of the upper theoretical bound on tti/jO 
(777/jO < ISOGeV) and the lower experimental bounds on the SUSY particles, the above 
condition cannot be fulfilled for any pair of SUSY particles (mi -|- 7772) and these loops 
will not contribute either. Thus, the only radiative corrections that will contribute are the 
self-energy loops on the iJ+-line and all vertex corrections. 

According to the particles in the loops, we have radiative corrections with sneutrinos 
and staus (Fig. la with sleptons), with stops and sbottoms (Fig. la with squarks) and with 
charginos and neutralinos (Fig. lb). (Note that the radiative corrections with ordinary 
quarks and with Higgs bosons are CP invariant.) The radiative corrections with scalar 
quarks are proportional to 777^ and, in general, one would expect that they should give 
the main contribution to Swh°- This contribution is enhanced, in addition, by the colour 
factor 3 that multiplies each diagram with squarks. Having at our disposal the parameters 
of the scalar mass matrices, one can achieve masses of the squarks that are low enough 
to be kinematically allowed in the considered range of 777 777^+ > {m^^^ + 777^^)^, still 
respecting the experimental bounds. However, at such small masses of the squarks, the 
CP-invariant radiative corrections to the denominator in (1.1) also grow and one can no 
longer expect that our first order formula (2,10) would be a good approximation. The 
performed numerical analysis confirmed these arguments. That's why we shall not consider 
the contribution of loops with squarks. Physically, this means that we assume that they 
are heavy and the decay tb is not allowed kinematically. In the commonly discussed 
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models of SUSY breaking, the squarks are much heavier than sleptons, charginos and 
neutralinos. 

Thus, finally we are left with the loop corrections on Fig. la with sleptons and on Fig. 
lb with chargino and neutralinos. The full analytic expressions for SY^. are rather lengthy 
but they considerably symplify for ^e{6Y^^). 

For the loops with staus and sneutrinos we have: 

XedYfKfVW, = "iZT^t^^ ^ 9m(toI)„WL)9".(Ci"+cl"+C<") 

m=l,2 

C^') = Cx{m%,ml,,mlmlr4^,ml) (2.11) 

m,n=l,2 

CP = Cx{mjj+,m^w,mlml,mlmlj (2.12) 



^ei5YfP)iTV) = / , sin(/3-a) ^ {mj^ - mD^mUglUTZfj^m 



SS'^ = 5o(m|+,mi,mi) (2.13) 



The loops with charginos and neutralinos give: 



27r . _ 

i=l,2 
/c,i=l,2,3,4 



Cx^ = C'x(m^+,mff,,m|,m|o,m~+,m|o). (2.14) 



5Re(5y5^^)(x°X+X+) = £ E {^"^(/iiO,^, A, + fM,^i)\m^^^n{2C^^ + + 



fc=l,2,3,4 



+ m|^+9m(Cf + "ih^^"i(<^2 )] 
+ m~^m^-,^m{fg0^jA*^ + /^0^^.4)9m(c(') + ^ + ) 

+ m^om~+9m(/^,Oi^,I; + /^0|;.4)9m(ci') + ^ + )) 

+ m^om~+S>m(/^,Oi^,.4- + fgol^A*,)^m{cf^ + )} 

<^J^ = Cx(m|+,m^,m|,m^+,m|o,m2 + ). (2.15) 

A-j A-fc A.j 
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Figure 1: The 1-loop diagrams in MSSM with complex parameters that contribute to 5^^° 

k=l,2,3,i 

i?f = Bo{m]j+ , m|o , ml+ ) (2.16) 

Here the imaginary parts of the PV integrals Cx{rn?jj+ , ^n^, m|, mg, mf , ml) and 
So(m^+, TTig, mf ) enter. The appropriate analytic expressions are given in Appendix B. 



3. Numerical Results 

Here we present our numerical analysis for the dependence of Swh° o^i the MSSM parame- 
ters. Taking into account the lower experimental and the upper theoretical bounds on the 
mass m^o, we consider the range 96 < rrij^o < 130 GeV. Our analysis showed a very weak 
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Table 1: The masses of the superparticles for the parameters (3.1) and (pr = (j^i = 



dependence on m^jO and the results here are presented for 

ruho = 125 GeV. 

As explained in Section 2 the diagrams with squarks cannot be considered within the 
a^-approximation used here and our analysis will be based on the diagrams on Fig. la 
with sleptons, and on Fig. lb. Thus our numerical results will depend on the parameters 
of the slepton, chargino and neutralino sectors. In order not to vary too many parameters 
we fix part of the SUSY parameter space: 

M2 = 250 GeV, Me = Ml -5 GeV, Ml = 120 GeV, 

\Ar\ = 500 GeV, = 150 GeV. (3.1) 

Assuming the GUT relation only for the absolute values of Mi and M2 (|Mi| = 
I tan 9w 1-^2 1), we keep ^1, the phase of Mi, as a physical phase. The phase of M2 can be 
rotated away. The other phase relevant for our considerations is the phase (pj- of A.,-- The 
varied parameters thus, are the charged Higgs mass, tan/3 and the CP-violating phases. 
We consider tan /? in the interval 

3 < tan/3 < 9. 

As it is well known, in order Siyi^o to be nonzero we need both new decay channels 
opened and CP violating phases. In accordance with this we have three cases: i) When 



tan(/3)=9 
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Figure 2: Swh" a-s a function of a) for different values of tan/3, Ml = 120 GcV; solid lines 
are for 0^ = — 7''/2, ^1 — 0; dashed lines are for 0^ = 0, (/)i = —tt/2. b) for different values of Ml, 
at tan/3 = 9, (/>r = —tt/2, (f)i = 0. 
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phases 



(a) (b) 

Figure 3: Syirho at tan/? = 9 versus the CP violating phases (f)r and (f)! for different values of 

<Pn [ (pii ^ 0, 7r/4, 7r/2]. The solid lines are for (j)^ = [~7r,7r] while <j>i = 0; the dashed lines are for 
01 = [— TT, tt] while 0T = 0, a) for m//+ = 237 GeV (= mp + 'm.+ ) and b) for m//+ = 387 GeV 
(=m-++m^o). 

only the decay channels vf^ are open (Fig. la). Then the phase (pr is responsible 

for CP violation, ii) When the decay channels —>■ xtXk open only (Fig. lb). In this 
case CP violation is due to the phase , and iii) When both — i>f^ and H'^ — x/"x^ 
decay channels are kinematically allowed ( all diagrams of Fig.l). In this case the two 
phases (f)r and cjyi contribute. 

Examples of the asymmetry as function of for cases i) and ii) are shown on Fig. 2a 
for different values of tan/3. The relevant SUSY mass spectrum is presented in Table 1. 
Numerically the two cases are obtained taking cpr 7^ 0, ^i=0 for case i), and 0r = 0, / 
for case ii)^. It is clearly seen that in both cases the asymmetry strongly increases with 
tan (3. The positions of the spikes correspond to the threasholds of the decay modes to the 
intermediate particles in the loops. It is seen that for mjj+ < 300 GeV the asymmetry is 
dominated by the light staus and sncutrinos in the SUSY mass spectrum - the solid lines 
on Fig. 2a, while for 111^+ > 300 GeV the asymmetry is determined by the charginos and 
neutralinos - the dashed lines on Fig. 2a. In both cases, i) and ii), the asymmetry reaches 
up to 10~^. The dependence on Ml for case i) is seen on Fig.2b. Case iii), when all relevant 
SUSY particles can be light, is described by the algebraic sum of the two graphs at a given 
tan /3 and we don't present it separately. In all these cases the asymmetry does not exceed 
few percents. 

Up to now, all presented results are for the phase of ^ zero, c/)^ = 0, i.e. for positive 
//'s. Negative roughly speaking just flips the sign of (^^y/jO. The effect of a non zero phase 

^The two cases i) and ii) can, surely, be obtained alternatively varying the mass parameters M2, /i and 
Ml. 
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4>^ is seen on Fig. 3, where we show the dependence of the asymmetry on (j)r (solid lines), 
and on 0i (dashed lines) for = 0, 7r/4 and 7r/2. The values of are chosen to be 

near the threasholds of the decay channels: mij+ = mc, + m-+ = 237 GeV on Fig.3a, and 

mfj+ = rn-+ + m^o = 378 GeV on Fig.3b. In all cases a CP violating phase of u does 

X2 xi 

not change the form of the curves but rather shifts the positions of the maximum and, 
in general, increases the absolute value of the asymmetry. Note that even for = the 
maximal effect is not achieved for 0r(0i) = — 7r/2. 

4. Conclusions 

We have considered the CP violating asymmetry 5\Yho of the decay rate difference between 
W^hP and W^h^ induced by one loop radiative corrections in MSSM with 

complex parameters. This decay is important for relatively low and tan /?. This in 

turn determines the importance of Siyj^o only if there are relatively low SUSY masses. We 
have considered the contribution from mn and mf, and/or and m^o in the loops, and 
thus the sensitivity to the phases of At, Mi and /j,. We work in first order of the weak 
coupling constant a«, and this approximation is not enough to consider the contribution 
of light stops and sbottoms. Typical values for the asymmetry are about 10~2 ^ 10-3, 
the main contributions being from D and f for mfj+ < 300 GeV, and from x'^ and for 
™-_ff+ > 300 GeV. The dependence on different values of tan P is examined. 

The approximate number of H'^^s needed to measure (5^^^° is Nff± > l/[(5^^o Br{H'^ 
W^hP)], which, for ^^rf^o ~ 10~^ and a branching ration 2± 10%, implies Nfj± > 10^. 

Charged will be produced at the Tevatron in FermiLab if mfj+ < 300 GeV, and 
at LHC in CERN if mfj+ < 1000 GeV. As the cross sections for production at pp and 
pp collisions decreases strongly for low tan/3, the required number Nfj± is too large for 
the planned luminosities at the hadron colliders [9] . For example, at LHC with integrated 
luminosity per year L = 100 fb~^, at mH+ = 500 GeV and tan/3 = 10, with an efficiency 
for the signal e = 2,6%, using the results of [9], for the ratio of the signal (S) and the 
background {B) events, we obtain S/B 0:^ 65/3770. 

More promising are the linear e^e~ colliders. In this case the charged Higgs will be 
copiously produced, the main production mechanism being e+e^ H^H~. Thus, the 
only parameter for the production cross section, at tree level, is the Higgs boson mass 
mH+. For a collider at -^/s = 800 GeV with luminosity L = 500 fb"^, the cross section 
is ~ 29 fb for mH+ = 200 GeV, and ~ 12 fb for mH+ = 300 GeV, which corresponds to 
1.5 X 10^ and 6 x 10^ H^H~-pairs, respectively [10]. For the CLIC collider, at ^/s = 3 
TeV, the cross section is 3 fb for mfj+ = 400 GeV which, for L = 800 fb~^, corresponds 
to 2.4 X 10^ charged Higgs pairs. This implies that at the NLC higher luminosities will be 
needed for such an asymmetry to be measured. 
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A. Lagrangian with complex couplings 



The mass matrices and their diagonahzation matrices are defined in the Appendicies A of 
[1] - for charginos and neutrahnos and of [2] - for staus and sneutrinos. Here we give only 
the pieces of the interaction Lagrangian we use. 

1. Lagrangian with neutral Higgses 

m,n=l,2 

jC^hOuu = 9 ^ sin(a + /3) uIVl (A.2) 

/ cos 

^ho^fxl = 9 Xi {AU Pl + Pr) xl (A.3) 
^hO^Ut = 3 xt {a*j Pl + Aji Pr) (A.4) 
^howw = gmw sm{P - a)W+W-'^h^ (A.5) 

where 

Cmn = (iLT^Lm^Ln + (^RR^RmT^Rn + (^RLT^RrrJ^Ln + (^R'^LrrJ^Rn 

ruz , 1 . 9 /, x . / ^^ ml sin a 

Cll = ^(- o + ew) sm a + /3 + -, 

cos Uw 2 mw cos p 

Crr = sm Ow sm(a + /3) H -, 

cos Uw fnw cos p 

c%T, = — - (u cos a + A*, sin a) 

^RL = o ^L. a + Ar sin a) = cI*r. 
zmw cos p 

Aik = ^ (sin a Qi^ + cos a S'l'^), 

Q'lk = \[Ni3{Nk2 - Nki tangly) + {I < — > k)], 

S'l'k = l[Ni4{Nk2 - Nkitanew) + {I ^ k)], 

Aij = sin a Qij — cos a Sij, 

Qij = --j=Ui2Vji, Sij = -^UaVj2. 

2. Lagrangian with charged Higgses^ 

^h'^H+w- = -f cos(a - P) {{H- da - {H+ da ^-'"} (A.6) 

^Thc correspondence with the notation in refs.[l] and [2] is gf^''^ = F^'^, g/{\/2mw)g4 = Gl, 

g/iV^mw)gl = Gi 
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^H+xtxl = 9 {H-Xk ( fk;PL + ft; Pr ) Xt + H+^t ( fjkPR + fik Pl ) xt] (A.7) 



where 



fR 



kj 



sin P 
COS /3 



NmVji + ^(^fe2 + Nki tanew)Vj2 



(A.8) 



^LL — ""^^ tan f3 — myy sin 2/3 
^LR = "^r (^T tan 13 + ji) 

3. Lagrangian with 

(A.9) 
(A.IO) 



a 



Here 



B. Absorptive parts of the Integrals 

If we use the notation: 

P° = g2 _ ml = {q+ pjf - ml 

then the Passarino-Veltman two- and three-point functions [11], when the loop integrals 
are given in 4-dimentions, are: 

1 r 1 

(B.l) 

(B.2) 
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ZTT^ j 
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ivr^ J 





For our diagrams in the considered iJ"^ decay we have: 

Pi=Ph+ ^Pi='m]j+, p2=ph^pl = ml, ipi-P2f = m^r- 
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Then for the the absorptive parts of the integrals, when the particles with masses mo and 
nil are put on mass shell we obtain: 



27r| k 



Here 



I k I 

X{x,y,z) 
A 
a 



SJm Co{m'jj+,mw, ml, ml, ml, ml) 

r< I 2 2 2 22 2\ 

ksmCi{mjj+,myir,mi^,mQ,mi,m2) 
Qm C2 {m]j+ , my^ , ml, ml, ml, ml) 



■ n 



In I 



a + b, 



"^rriH+lphl a-b 
ml A- ipH+Ph) B 
A 

-{PH+Ph)A + m]j+B 
A 



A = 



(PH+Ph) = 



^H+ml - (PH+Ph)'^, 

2mH+ 

+ + — 2xy — 2xz — 2yz 



mjj+ + mf^ 



m 



w 



\Ph\ 



\^''^{m]j+,ml,m]^) 
2mH+ 



nk^ ^ ,a + 6, 
m 



B 



2\Ph\ '"'a-b ' 2mH+ \Ph 

ml + ml-ml + 2plk\ b=-2\t\\rh\, 



k'pl-l 



In 



a + b, 



(B.4) 
(B.5) 
(B.6) 
(B.7) 



2\k\\rh\ 



mT — m. 







m 



H+ 



2mH+ 



pI 



mjj+ -m^ + ml 
2mH+ 
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